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ABSTRACT 

Context. We discuss here the properties of the double cluster Abell 1758, at a redshift z~0.279, which shows strong evidence for merging. 
Aims. We analyse the optical properties of the North and South clusters of Abell 1758 based on deep imaging obtained with the CFHT archive 
Megaprime/Megacam camera in the g' and r' bands, covering a total region of about 1.05x1.16 deg'^, or 16.1 x 17.6 Mpc*. Our X-ray analysis is 
based on archive XMM-Newton images. Numerical simulations were performed using an N-body algorithm to treat the dark matter component, 
a semi-analytical galaxy formation model for the evolution of the galaxies and a grid-based hydrodynamic code with a PPM scheme for the 
dynamics of the intra-cluster medium. We have computed galaxy luminosity functions (GLFs) and 2D temperature and metallicity maps of the 
X-ray gas, which we then compared to the results of our numerical simulations. 

Methods. The GLFs of Abell 1758 North are well fit by Schechter functions in the g' and r' bands, but with a small excess of bright galaxies, 
particularly in the r' band; their faint end slopes are similar in both bands. On the contrary, the GLFs of Abell 1758 South are not well fit 
by Schechter functions: excesses of bright galaxies are seen in both bands; the faint end of the GLF is not very well defined in g'. The GLF 
computed from our numerical simulations assuming a halo mass-luminosity relation agrees with those derived from the observations. From 
the X-ray analysis, the most striking features are structures in the metal distribution. We found two elongated regions of high metallicity in 
Abell 1758 North with two peaks towards the center. On the other hand, Abell 1758 South shows a deficit of metals in its central regions. 
Comparing observational results to those derived from numerical simulations, we could mimic the most prominent features present in the 
metallicity map and propose an explanation for the dynamical history of the cluster. We found in particular that in the metal rich elongated 
regions of the North cluster, winds had been more efficient in transporting metal enriched gas to the outskirts than ram pressure stripping. 
Results. We confirm the merging structure of each of the North and South clusters, both at optical and X-ray wavelengths. 
Conclusions. 

Key words. Galaxies: clusters: individual (Abell 1758), Galaxies: luminosity function 



X 



1. Introduction 

Environmental effects are known to have an influence on 
galaxy evolution, and can therefore modify Galaxy Luminosity 



* Based on archive data retrieved from the Canadian 
Astronomy Data Centre Megapipe archive and obtained with 
MegaPrime/MegaCam, a joint project of CFHT and CEA/DAPNIA, 
at the Canada-France-Hawaii Telescope (CFHT) which is operated by 
the National Research Council (NRC) of Canada, the Institut National 
des Sciences de I'Univers of the Centre National de la Recherche 
Scientifique (CNRS) of France, and the University of Hawaii. The 
X-ray analysis is based on XMM-Newton archive data. This research 
has made use of the NASA/IPAC Extragalactic Database (NED) 
which is operated by the Jet Propulsion Laboratory, California 
Institute of Technology, under contract with the National Aeronautics 
and Space Administration, and of the SIMBAD database, operated at 
CDS, Strasbourg, France. 



Functions (hereafter GLFs). This is particularly obvious in 
merging clusters, where GLFs may differ from those in non- 
merging (relaxed) clusters, and where GLFs may also be ob- 
served to differ be tween one photometric band and another (e.g. 



Boue et al. Il2008h and references therein). GLFs also allow to 



trace the clust er formation history, as shown for example in the 
case of Coma jAdami et al. I2OO7I) . 

This dynamical history can also be derived by analyzing 
the temperature and metallicity distributions of the X-ray gas 
in clusters. Such maps have revealed that in many cases, clus- 
ters with emissivity maps showing a rather relaxed appearance 
could have very disturbed temperature and m etallicity distri- 



butions (see e.g. iDurret & Lima Neto I l2008l and references 
therein), meaning that they have undergone one or several 
mergers in the last few Gyrs. The study of the thermal structure 
of the intra-cluster medium (ICM) indeed provides a very inter- 
esting record of the dynamical processes that clusters of galax- 
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ies have experienced during their formation and evolution. The 
temperature distribution of the ICM gives us insight into the 
process of galaxy cluster merging and on the dissipation of the 
merger energy in form of turbulent motion. Metallicity maps 
can indeed be regarded as a record of the integral yield of all 
the different stars that have released their metals through super- 
nova explosions or winds during cluster evolution. 

The comparison of temperature and metallicity maps 
to the results of hydrodynamical numerical simulations al- 
low to characterize the last merging events which have 
taken or are taking place. For example, the comparison of 
the complex temperature and meta llicity maps of Abell 85 
Durret. Lima Neto & Forman (120051) with the numerical simu- 
lations bv lBourdin et al. 1(120041) show that two or three mergers 
have taken place at various epochs in this cluster in the last few 
Gyrs, besides the ongoing mer ger seen as a filame nt made of 



groups falling on to the cluster (IDurret et al. I2OO3I 



We have become interested in pairs of clusters, where 
the effects of merging are expected to be even stronger In 
some cases, one of the clusters shows itself a double structure 
(Abell 223, Abell 1758 North). By coupHng deep optical multi- 
band imaging with X-ray r naps, we have recen tly analyzed the 



Abell 222/223 cluster pair (iDurret et al. 112010b . We found that 



Abell 222 (the less perturbed and less massive cluster) had 
GLFs well fit by a Schechter function, with a steeper faint end 
in the r' band than in the g' band, implying little star forma- 
tion; its X-ray gas showed quite homogeneous temperature and 
metallicity maps, but with no cool core, suggesting that some 
kind of merger must have taken place to suppress the cool core. 
This was confirmed by the distribution of bright galaxies in 
this cluster, which also suggests that this cluster is not fully re- 
laxed. Abell 223 (the most perturbed and massive cluster), was 
found to have comparable GLFs in both bands, with an excess 
of galaxies over a Schechter function at bright magnitudes. Its 
temperature and metallicity distributions were found to be very 
inhomogeneous, implying that it has most probably just been 
crossed by a smaller cluster, which now appears at the north 
east tip of the maps. Note th at a bridge of galaxies seems to ex - 
ist between the two clusters dDietrich. Clowe & Soucail l2002h . 
as well as a possible dark matter filament joining the two clus- 
ters bietrichet aTl (|2005|) . 



The Abell 1758 cluster, at a redshift of 0.279, was anal- 



ysed in X-rays bv lDavid & Kempnerl ( 120041) based on Chandra 
and XMM-A^evvfo« data. These authors showed that this cluster 
is in fact double, with a North and a South component sepa- 
rated by approximately 8 arcmin (2 Mpc in projection), both 
undergoing major mergers, with evidence for X-ray emission 
between the two clusters. However, very little has been pub- 
lished on this system at optical wavelengths, and few galaxy 
redshifts are available. As a second study of cluster pairs, we 
chose to analyse this system by coupling archive optical CFHT 
Megacam data with archive XMM-Newton data. This allowed 
us to compute GLFs in two bands as well as temperature and 
metallicy maps for the ICM, which were compared to the re- 
sults of numerical simulations. 



F or a redshift o f 0.279, Ned Wright's cosmology calcula- 
toJil (IWrightll2006l) gives a luminosity distance of 1428 Mpc 
and a spatial scale of 4.233 kpc/arcsec, giving a distance mod- 
ulus of 40.77 (assuming a flat ACDM cosmology with Hq = 
70 km S-' Mpc-', Q.m = 0.3 and Qa = 0.7). 

The paper is organised as follows. We describe our optical 
analysis in Section 2, and results concerning the observed and 
simulated galaxy luminosity functions in Section 3. The X-ray 
data analysis and results, including temperature and metallic- 
ity maps, are presented in Section 4. The results of numerical 
simulations that were run to help us to account for the X-ray 
temperature and metallicity maps are described in Section 5. 
An overall picture of this cluster pair is drawn in Section 6. 



2. Optical data and analysis 

2.1. The optical data 

We 



(Gwvn 



hav e retr ieved 
2009h 



from the CADC Megapipe archive 
the reduced and stacked images in the 
g' and r' bands (namely G008.203.1404-50.518.G.fits and 
G008.203.1404-50.518.R.fits) and give a few details on the ob- 
servations in Table [1] Observations were made at the CFHT 
with the Megaprime/Megacam camera, which has a pixel size 
of0.186x 0.186 arcsec-. 

Table 1. Summary of the observations. 



Filter 


g' 


r' 


Number of coadded images 


4 


9 


Total exposure time (s) 


1800 


4860 


Seeing (arcsec) 


0.75 


0.65 


Limiting magnitude (5a) 


27.1 


26.8 



We did not use the catalogues available for these images, 
because they were made without masking the surroundings of 
bright stars, so we preferred t o build masks first, then to ex- 
tract sources with SExtractor (Bertin & Arnoutsi 1996 1. The 



total area covered by the images was 20403x22406 pixels^, or 
1.05x1.16 deg- (16.1 x 17.6 Mpc^ at the cluster redshift). 

Objects were detected and measured in the full r' image, 
then measured in the g' image in double image mode (i.e. the 
objects detected in r' were then measured in g' exactly in the 
same way as in r'). Magnitudes are in the AB system. The ob- 
jects located in the masked regions were then taken out of the 
catalogue, leading to a final catalogue of 298,170 objects. We 
created masks around bright stars and image defects in the por- 
tion of the image covered by Abell 1758 and by the ring around 
the cluster that was used to estimate the background galaxy 
contamination to the GLF (due to its relatively high redshift 
the cluster does not cover the entire image). After taking out 
masked objects we were left with a catalogue of 286,505 ob- 
jects with measured r' magnitudes, out of which 277,646 also 
have measured g' magnitudes. 



http://nedwww.ipac.caltech.edu/ 
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Fig. 1. Central surface brightness in the r' band as a function 
of r' magnitude. The horizontal and oblique full lines isolate 
the star sequence (below the horizontal line and between the 
oblique lines). The two vertical lines correspond to r' - 17.75 
where the stars stop being saturated, and r' - 22, where the fit 
to calculate the star sequence was limited (see text). 

Since the seeing was better in the r' band (see Table [T]), we 
performed our star-galaxy separation in this band. 

2.2. Star-galaxy separation 

In order to separate stars from galaxies, we plotted the maxi- 
mum surface brightness /i„„„ in the r' band as a function of r'. 
The result is shown in Fig.[T] 

The best fit to the star sequence visible on Fig.[T]calculated 
for 17.75 < r' < 22 is ju„,„ = 0.996r' - 0.524, with standard 
deviations on the slope and constant of 0.002 and 0.035 re- 
spectively. The point-source (hereafter called "star") sequence 
is clearly visible for r' < 22, with the star saturation showing 
well for r' < 17.75. We will define galaxies as the objects with 
fimaAf') > 17.4 for r' < 17.75, and as the objects above the line 
of equation ^,„„. = 0.996/ - 0.274 for r' > 17.75. Stars will 
be defined as all the other objects (see Fig.[TJ. The small cloud 
of points observed in Fig. [1] under the star sequence is in fact 
defects, but represents less than 2% of the number of stars. We 
thus obtained a star and a galaxy catalogue. 

As a check to see up to what magnitude we could trust 
our star-galaxy separation, we retrieved the star catalogue from 
the Besanfon model for our Galaxy (Robin et al. 2003) in a 
1 deg- region centered on the position of the image analysed 
here. Such a catalogue is in AB magnitudes (as ours) and is 
corrected for extinction. In order for it to be directly compara- 
ble to our star catalogue, we corrected our star catalogue (and 
our galaxy catalogue as well, for later purposes) for extinction: 
0.0531 mag in g' and 0.0385 mag in r' (as derived from the 
Schlegel et al. 1998 maps). 




RA (deg) 



Fig. 2. Positions of the objects (stars and galaxies) in the 
Abell 1758 North (red) and South (green) catalogues. The ob- 
jects used to estimate background counts (see text) are shown 
in magenta. Note that the figure covers 1x1 deg^, and is slightly 
smaller than the images. 



The r' magnitude histogram of the objects classified as stars 
in our r' image roughly agrees with the Besan§on star catalogue 
for r' <22. However, for r' >21, we start to detect more stars 
than predicted by the Besangon model (the difference is only 
about 15% at r' - 21.5 but becomes 25% in the r' = 22 bin, 
and the difference continues to increase at fainter magnitudes). 

We will therefore consider that our star-galaxy separation 
is correct for r' < 22. For fainter magnitudes, we will compute 
galaxy counts by counting the total number of objects (galaxies 
plus stars) per bin of 0.5 mag, and considering that the number 
of galaxies is equal to the total number of objects minus the 
number of stars predicted in each bin by the Besanfon model. 

We extracted from the star and galaxy catalogues two cat- 
alogues as large as possible corresponding to the North and 
South clusters. Their respective positions were taken to be 
the X-ray centers of the two clusters, as derived from XMM- 
Newton data: 203. 1851, +50.5445 (J2000.0, in degrees) for the 
North cluster, and 203.1335,+50.4138 for the South one. The 
maximum possible radius to obtain independent catalogues for 
each of the two clusters was 0.0675 deg, or 1.03 Mpc at a red- 
shift of 0.279. 

We obtained for each of the two clusters three complemen- 
tary catalogues (with g' and / magnitudes): objects classified 
as galaxies (classification valid at least for r' < 22), objects 
classified as stars, and a complete catalogue of galaxies+stars 
which will be used for r' > 22. The positions of the galaxies in 
the regions of the two clusters are shown in Fig. |2] 
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Fig. 3. Point source completeness as a function of magnitude in 
percentages for Abell 1758 North (top) and South (bottom) in 
g' (left) and r' (right) for point-like objects (see text). 



2.3. Catalogue completeness 

The completeness of the catalogue is estimated by simulations. 
For this, we add "artificial stars" (i.e. 2-dimensional Gaussian 
profiles with the same Full-Width at Half Maximum as the av- 
erage image Point Spread Function) of different magnitudes to 
the CCD images and then attempt to recover them by running 
SExtractor again with the same parameters used for object de- 
tection and classification on the original images. In this way, 
the completeness is measured on the original images. 

In practice, we extract from the full field of view two 
subimages, each 1300 x 1300 pixels-, corresponding to the po- 
sitions of the two clusters on the image. 

In each subfield, and for each 0.5 magnitude bin between 
r' = 20 and 27, we generate and add to the image one star 
that we then try to detect with SExtractor, assuming the same 
parameters as previously. This process is repeated 100 times 
for each of the two fields and bands. 

Such simulations give a completeness percentage for stars. 
This is obviously an upper limit for the completeness level for 
galaxies, since stars are easier to detect than galaxies. However, 
we have shown in a previous paper that this method gives a 
good estimate of the complet eness for normal ga laxies if we 
apply a shift of ~ 0.5 mag (see lAdami et al. l2006h . Results are 
shown in Fig. [3] 

From these simulations, and taking into account the fact 
that results are worse by ~ 0.5 mag for mean galaxy popula- 
tions than for stars, we can consider that our galaxy catalogues 
are complete to better than 80% for g' < 25.8 and r' < 25.6 in 
both clusters. 



2.4. Galaxy counts 

The surfaces covered by the Abell 1758 North and South cata- 
logues (after excluding masked regions) are 0.01412 deg^ and 



0.01418 deg^ respectively. Galaxy counts were computed in 
bins of 0.5 mag normalized to a surface of 1 deg^. 

For r' < 22, galaxy counts were derived directly by com- 
puting histograms of the numbers of galaxies in the Abell 1758 
North and South catalogues. For r' > 22, we built for each clus- 
ter histograms of the total numbers of objects (galaxies +stars) 
and obtained galaxy counts by subtracting the numbers of stars 
predicted by the Besanfon model. The resulting galaxy counts 
will be used in the next section to derive the GLFs for both 
clusters in both bands. 

As a test, we considered the galaxy counts in the 21.5-22.0 
magnitude bin computed for both clusters with the two meth- 
ods (i.e. first method: considering that the star-galaxy separa- 
tion is valid, and second method: considering the total num- 
ber of objects (galaxies+stars) and subtracting the number of 
stars predicted by the Besangon model to obtain the number of 
galaxies). In all cases, the differences are smaller than 4%. 

Note that no k-correction was applied to the galaxy magni- 
tudes. 

3. Results: colour-magnitude diagrams and galaxy 
luminosity functions 

In order to compute the galaxy luminosity functions of the 
two clusters, we need to subtract to the total galaxy counts the 
number counts corresponding to the contamination by the fore- 
ground and background galaxies. 

For galaxies brighter than r' - 22 we will select galax- 
ies with a high probability to belong to the clusters by draw- 
ing colour-magnitude diagrams and selecting galaxies located 
close to this relation. A few spirals may be missed in this way, 
but their number in any case is expected to be small, as ex- 
plained in Sect. 3.1 (also see e.g. Adami et al. 1998). For galax- 
ies fainter than r/ = 22 we will subtract galaxy counts statisti- 
cally. 

3.1. Colour-magnitude diagram 

The g' - r' vs. r' colour-magnitude diagram is shown in Fig.|4] 
for the two clusters together (the diagrams are the same for 
both clusters). A sequence is well defined for galaxies in the 
magnitude range 17 < r' < 22 in both clusters. We computed 
the best fit to the g' - r' vs. r' relations in this magnitude range 
by applying a linear regression. We then eliminated the galaxies 
located more than 3cr away from this relation and computed the 
g' - r' vs. r' relation again. 

The equation of the colour-magnitude relation is found to 
be: g' - r' - -0.0436r' + 2.108 with an r.m.s. on the constant 
o- = 0.11. 

We also plot in Fig.|4]galaxies with measured spectroscopic 
redshifts. We can see that the positions of the galaxies belong- 
ing to the cluster according to their spectroscopic redshifts, i.e. 
with redshifts in the [0.264, 0.294] interval fall very close to 
the best fit to the colour-magnitude relation. 

In view of this, for r' < 22, we will consider hereafter that 
all the galaxies located within +1.5cr of the colour-magnitude 
relation (i.e. between the two black lines of Fig. |4]i belong to 
the cluster. 
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Fig. 4. ig'-rf) vs. r' colour-magnitude diagram for Abell 1758 
(North and South) for objects classified as galaxies from the 
P^max - r' magnitude relation. The vertical dashed lines indicate 
the magnitude interval where the colour-magnitude relation 
was computed. The long oblique dashed line shows the mean 
colour-magnitude relation; the short oblique dotted Unes indi- 
cate intervals of +lcr, 2cr and 3cr around the colour-magnitude 
relation. The interval finally adopted of 1 .5cr is shown with full 
lines. The filled circles show the galaxies with measured red- 
shifts, colour-coded as follows: black: galaxies belonging to 
the cluster according to their spectroscopic redshifts, cyan and 
magenta for galaxies with smaller and larger redshifts than the 
cluster, but inside the two circles where the two clusters were 
extracted, blue and red: galaxies with smaller and larger red- 
shifts than the cluster, but outside these two circles. 



We can note that the scatter cr = 0.11 is somewhat larger 
than found in the Abell 222/223 clusters at z-Q.2\ (Durret et al. 
2010), but the interval chosen for cluster membership (+1.5cr) 
remains smaller than that used for example by Lagana et al. 
(2010) in the redshift interval [0. 1 1,0.23]. 

The initial sample of galaxies comprises 1599 and 1438 
galaxies in the North and South clusters respectively (with no 
magnitude limit). Within the +1.5cr interval along the red se- 
quence, we are left with respective numbers of galaxies of 477 
and 340. These numbers become 192 and 116 galaxies for 
r' < 22. The North cluster is obviously richer than the South 
one. 

With this rather strict criterium, we obviously select galax- 
ies with a high probability to belong to the cluster, but we may 
lose some galaxies, in particular blue cluster galaxies falling 
under the sequence. 

We have therefore estimated the number of blue cluster 
galaxies lost by selecting galaxies within +1.5cr of the red se- 
quence in the following way. First, we computed histograms 
of numbers of galaxies within +1.5cr of the red sequence and 
below this sequence in bins of 1 absolute magnitude in the r' 
band. These counts were made in absolute magnitude bins to 
be comparable to the counts estimated from luminosity func- 
tions of field galaxies. The bins of interest here are between 
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Fig. 5. Galaxy counts in the g' (left) and r' (right) bands for 
magnitudes r' > 22 where the background must be subtracted 
statistically, in logarithmic scale. The counts in Abell 1758 
North and South are drawn in red and green respectively. The 
magenta dashed lines show the galaxy counts from the "local" 
background extracted in the annulus shown in Fig. |2] Error bars 
are Poissonian and are not plotted for clarity. 



Mr' = -22 and -19, roughly corresponding to r' - 18.5 and 
21.5. We then computed the number of foreground galaxies ex- 
pected. Since the comoving volume at z-Q.219 is 5.834 Gpc^ 
(Wright 2006), and each of our clusters covers an area of 
0.01431 deg^ on the sky, the volume in the direction of each 
cluster is 2024 Mpc^. By using the R band luminosity func- 
tion by Ilbert et al. (2005) in the 0.05-0.20 redshift bin (see 
their Fig. 6 and Table 1), we find that the percentages of "lost" 
galaxies are of the order of 30% for M,- = -19, of 10%-15% 
for M,-' = -20 and -21, and less than 10% for brighter galax- 
ies. 



3.2. Comparison field 

In order to perform a statistical subtraction of the background 
contribution for r' > 22, and since the cluster is quite distant 
and does not cover the whole field, we extracted background 
counts in an annulus surrounding the clusters (see Fig. |2). The 
annulus was centered on the middle position between the two 
clusters (203.1593,-h50.4792 J2000.0 in degrees), with an in- 
ner radius of 0.3232 deg (4.9 Mpc), and an outer radius of 
0.4444 deg (6.8 Mpc). The unmasked surface of the annulus 
is 0.278 degl 

We checked if there could be any contamination of 
background counts in this annulus by a group or clus- 
ter, and found a structure west-southwest of Abell 1758. 
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This appears as a cluster in SimbacH with coordinates 
13''34'««45.41', +50°26'01.4" (J2000.0) and redshift 0.085. 
The mean redshift for the 20 galaxies extracted from NED 
(http://nedwww.ipac.caltech.edu/) in this region is 0.0869, with 
a dispersion in redshift = 0.0015. This cluster falls just out- 
side the annulus which was used to estimate the background 
counts subtracted to the galaxy counts to compute the GLF, so 
its presence should not modify galaxy counts inside the annu- 
lus. 

We can note from Fig.|5]that galaxy counts are comparable 
in Abell 1758 North and South in the g' band, but differ in 
the r' band, where Abell 1758 North has more galaxies in the 
22 < r' < 24 magnitude range. 

3.3. Galaxy luminosity functions 

The Galaxy Luminosity Functions (GLFs) of Abell 1758 North 
and South were calculated in bins of 0.5 mag and normalized 
to 1 deg^. We subtracted the background contribution using as 
background galaxy counts the "local" counts in the annulus 
shown in magenta in Fig. |2] 

The GLFs are displayed for Abell 1758 North and South in 
Figs.|6]and|7]respectively (red points). The error bars drawn in 
these figures were taken to be 4 times the Poissonian errors on 
galaxy counts, as derived from detailed sim ulations previously 



A1758N 



20081 



performed by our team for similar data (see iBoue et al. 
Fig. 5). 

The GLFs (as a function of absolute magnitude) were fit by 
a Schechter function: 

S(M) = 0.4 In 10 (/.*/'+' e"-'" 
withy = 100-4 (M--M)^ 

The parameters of the Schechter function fits of the GLFs 
are given in Table|2] The absolute magnitude ranges considered 
are indicated for each fit, and GLFs and their fits are drawn in 
Figs.|6]and|2l 

If we look at Abell 1758 North (Fig. |6ll we see that a 
Schechter function fits rather well most of the GLF points for 
both bands. However, there is an excess of galaxies over a 
Schechter function in the very brightest magnitude bins, spe- 
cially in the r' band. There is also a "bump" around r' abso- 
lute magnitudes -17 and -17.5 which has no obvious explana- 
tion. Except for this feature, the GLFs are quite similar in both 
bands, and the faint end slopes are quite flat: a = -0.85. 

On the other hand, for Abell 1758 South (Fig.|2), the GLFs 
are not well fit by Schechter functions. In the g' band, there is 
a strong excess of galaxies at bright magnitudes, and the faint 
end slope is quite flat [a = -0.59), implying that star formation 
is weak in faint galaxies of the South cluster. In the r' band, 
the GLF also shows an excess at very bright magnitudes, and 
a strong dip for Mi- ~ -17.5. This dip is also seen in the g' 
band around -17, though it is not as pronounced as in the r' 
band. Note that in the South cluster the GLFs in the g' and r' 
bands have quite different faint end slopes: a - -0.59 in g' and 
a = -1.11 in r' . This lack of blue faint galaxies could suggest 
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Fig. 6. Galaxy luminosity functions for Abell 1758 North in the 
g' (top) and r' (bottom) bands, in logarithmic scale. The blue 
and red points correspond to the two galaxy selections (see 
text), and the best Schechter function fits are drawn with the 
same colours as the corresponding points. Note that at blight 
magnitudes, the points exactly coincide, and since the blue 
points were plotted after the red ones, they appear blue. Error 
bars are 4 times the Poissonian errors on galaxy counts (see 
text). 



that star formation has been quenched by a process linked to 
the merger, but could also be an artefact due to the method used 
here to derive the GLF (see end of Sect. 3.4). 

The fact that both clusters have GLFs diff'ering from simple 
Schechter functions is most probably due to the fact that both 
are undergoing merging p rocesses, as already pointed out by 



David & Kempner I (120041) from their X-ray study. We will dis- 



http://simbad.u-strasbg.fr/simbad/ 



cuss these results in the next Section when considering the tem- 
perature and metallicity distributions of the X-ray gas. These 
maps confirm that both clusters are indeed structures which are 
strongly pertubed by several mergers, so it is not surprising to 
see effects on the GLFs. 

Altogether, the GLFs in Abell 1758 do not strongly dif- 
fer from those derived in other clusters. The bright parts of 
the GLF Schechter fits (M,- < -19) for both clusters are quite 
similar in shape to th e GLFs recently obt ained by other au- 
thors (see for example I Andreon et al. 120081) . And the faint end 
slopes are within the broad range of values estimated by previ- 
ous authors for different clusters, cluster regions and photomet- 
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Table 2. Schechter parameters for galaxy luminosity functions. The first set of fits corresponds to galaxies selected from the 
colour-magnitude relation for r' < 22 and to a statistical background subtraction with galaxy counts taken from the CFHTLS 
Deep field counts for r' > 22. The second set of fits corresponds to galaxies selected from the colour-magnitude relation at all 
magnitudes. The last Une corresponds to the fit to the simulated GLF (see text). 



Cluster 


Filter 


Range 


O* 




M* 




a 




North 


8' 


-22.0,-16.5 


4384 ± 


462 


-20.06 ± 


0.12 


-0.86 ± 


0.05 




r' 


-23.5,-16.0 


4650 ± 


286 


-21.18 ± 


0.08 


-0.85 ± 


0.02 


South 


g' 


-22.0,-15.5 


5319 ± 


574 


-18.90 ± 


0.14 


-0.59 ± 


0.08 




r' 


-23.5,-16.0 


1582 ± 


229 


-21.34 ± 


0.16 


-1.11 ± 


0.04 


North 


s' 


[-22.0,-16.5] 


3244 ± 


236 


-20.33 ± 


0.09 


-1.00 ± 


0.02 




r' 


[-23.5,-16.0] 


3451 ± 


229 


-21.44 ± 


0.08 


-0.96 ± 


0.02 


South 


g' 


[-22.0,-15.5] 


1652 ± 


205 


-20.09 ± 


0.16 


-1.12± 


0.03 




r' 


[-23.5,-16.0] 


1849 ± 


195 


-21.18 ± 


0.14 


-1.06 ± 


0.03 


Simulated 




[-25.0,-15.0] 


755 ± 


21 


-22.09 ± 


0.04 


-1.01 ± 


0.01 



A175BS 




-20 

Magnitude i 



Al 758S 



10 1 

-24 




Fig. 7. Same as Fig.|6]for the South cluster. 



ric ba nds (see e.g. the compilation in Table Al of iBoue et al. 
20081) . except for the South cluster in the g' band, where the 
GLF seems unusually flat. 

Note that AbeU 1758 is at redshift 0.279, and few GLFs are 
available for clusters at such redshifts 



Andreon et al. 



(I2005h 

found more or less comparable faint end slopes of - 1 .03 and 
-1.30 for two clusters at redshifts ~ 0.3, but in the K band, so 
the comparison with our results is not straightforward. 



3.4. Galaxy luminosity functions only based on a 
colour-magnitude selection 

As a test to our method, we also derived the GLFs by consid- 
ering that all the galaxies located within +1.5cr of the colour- 
magnitude relation shown in Fig.|4]were cluster members, for 
all magnitudes. 

The corresponding points and GLF Schechter fits are 
shown in blue in Figs. |6] and |7] and the corresponding Schechter 
parameters are given in the second half of Table |2] We can see 
that for absolute magnitudes fainter than about -19 to -18 the 
data points of the GLFs start to differ The agreement between 
the Schechter fits based on the two methods is fair for the North 
cluster in both bands and for the South cluster in the r' band, 
though the Schechter parameters found in corresponding cases 
are not always within error bars. This strongly suggests that the 
errors on these parameters are underestimated, and this is prob- 
ably also the case for the errors on the GLF points themselves. 
On the other hand, the agreement is poor for the South cluster 
in the g' band. Therefore we cannot consider that the g' band 
GLF is well constrained in the South cluster 

These results clearly illustrate the difficulty to estimate 
GLFs, particularly at faint magnitudes, as already mentioned 
by a number of authors (see discussion in Durret et al. 2010). 
They also justify our choice not to attempt GLF fits with a 
higher number of free parameters, as would be obtained by fit- 
ting a gaussian at bright magnitudes plus a Schechter function 
at faint magnitudes, or two Schechter functions. 

3.5. Simulated galaxy luminosity function 

In connection with the numerical simulations presented in 
Section 5, we simulated a galaxy luminosity function using the 
halo mass-luminosity relation from Vale & Ostriker (2006). 
The halo mass for the simulated galaxies is taken from the 
galaxy formation model, which calculates the halo mass from 
the N-body simulation. All the simulated galaxies were put into 
25 magnitude bins between absolute magnitudes -25 and -15. 

The resulting GLF is shown in Fig. [8] We can see that it 
appears quite similar in shape to those shown in Figs. |6]and|7] 
A Schechter fit to this function is superimposed in Fig. |8]and 
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Fig. 8. Simulated luminosity function with best Schechter fit 
superimposed. 



its parameters are given in Table 2. We can see that although 
the value of M* is between half a magnitude and a magnitude 
brighter, the faint-end slope agrees well with the values derived 
from the observations. 



3.6. Do the North and South clusters have a cD 
galaxy? 

The images of the three brightest galaxies of each cluster are 
shown in Fig. |9l We can see that in the North cluster none of 
these galaxies is at the cluster centre, and there are obviously 
two subclusters. In the South cluster, none of the three brightest 
galaxies is perfectly in the centre either (the galaxy which is not 
circled is a foreground object). 

The similarity with the Abell 222/223 cluster pair (see 
Durret et al. II2OI0I figure 20) has led us to analyze the bright- 
ness profiles of the brightest galaxies in the Abell 1758 North 
and South clusters. We had found that Abell 223 (which re- 
sembles Abell 1758 North) had two Brightest Cluster Galaxies 
(BCGs), one of them showing a brightness profile decreasing 
slowlier than the other, and therefore resembling that of a cD, 
while the other one could be the central galaxy of an accreted 
group. 

The surface brightness profiles of the three brightest galax- 
ies in the two Abell 1758 clusters are displayed in Fig. [TO] 
We can see that the profiles of the two brightest galaxies of 
the North cluster decrease notably slowlier with radius that 
those of the other bright galaxies, and surprisingly, the pro- 
file of the second brightest galaxy is much flatter than that 
of the brightest one. This suggests that there are two domi- 
nant galaxies in the North cluster, quite similarly to what was 
seen in Abell 223, and this is another indication of a merger. 
The fact of having two dominant galaxies is probably another 
indi cation of a merger of two smaller systems, in agreement 
with lDavid & Kempnerl (120041) who suggested that at least two 
smaller clusters have crossed this North system. On the other 
hand, there is no dominant galaxy in the South cluster, where 
the profiles of the three brightest galaxies are comparable. 




Fig. 9. Optical images of Abell 1758 North (lop) and South 
(bottom). For each cluster, the brightest galaxy is circled in 
black, the second brightest in blue and the third brightest in 
cyan. The red circles show the cluster limits. The uncircled 
galaxy near the center of the South cluster is a foreground ob- 
ject. 



4. X-ray analysis 

We present in Fig.[TT]the X-ray surface brightness isocontours 
oveplotted on the optical r' band image of Abell 1758 and in- 
dicate the positions of the three brightest galaxies of each clus- 
ter. The three galaxies of each cluster are the remaining central 
galaxies from the ancient sub-clusters. Another important thing 
to notice is their position: the two brightest galaxies one near 
the other and the third one is radially opposite. This is an indi- 
cation of a recent merger where we clearly see that while the 
gas has already settled down, the brightest galaxies are not in 
the center of Abell 1758 North and A1758 South. We wiU re- 
visit this merging scenario in the following sections. 

4.1. Data reduction 

Abell 1758 was observed for -57 ks with XMM-Newton with 
the Medium filter inserted. The XMM-Newton ODF files were 
processed using SAS version v8.0. The MOS and pn files were 
filtered excluding all events with FLAG > and PATTERN 
> 12, and FLAG > and PATTERN > 4, respectively. Light 
curves were made in the 1-10 keV energy band and periods 
where the background value exceeded the mean value by more 
than 3cr were excluded. We considered events inside the field 
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Bright galaxy profiles 




50 100 150 200 

radius (px) 

Fig. 10. Surface brightness profiles of the three brightest galax- 
ies in both clusters (intensities are normalized to the central 
value and in arbitrary units, and radii are in pixels). Full and 
dashed lines correspond to galaxies in the North and South 
clusters respectively. For each cluster, the brightest galaxy is 
in black, the second brightest in blue and the third brightest in 
cyan. 



of view (FOV) and excluded all bad pixels. Flare filtering left 
live times of 16928 s, 12734 s and 16929 s in the MOSl, MOS2 
and pn cameras respectively. 

The background was taken into account by extracting 
MOSl, MOS2 and pn spectra from the pubUcly available EPIC 
blank sky templates of Andy Read (Read & Ponman 2003). 
The background was normalized using a spectrum obtained in 
an annulus (between 12.5-14 arcmin) where the cluster emis- 
sion is no longer detected. 

4.2. Spectrally measured 2D X-ray maps 

We performed quantitative studies using X-ray spectrally mea- 
sured 2D maps to derive global properties of these two clusters. 
These maps were made in a grid; for each spatial bin we set a 
minimum count number of 900 (after background subt raction) . 
For th e spectral fits, we used XSPEC version 11.0.1 (lArnaud 
1996h and modeled the obtained spectra with a MEKAL sin- 
gle tempe rature plasma emission model (bremsstrahlung + line 
emission 



Kaastra & Mewe 1993: Liedahl et al. 1995). The 



free parameters are the X-ray temperature (kT) and the metal 
abundance (metallicity). Spectral fits were made in the energy 
interval of 0.7-8.0 keV with the hydrogen column density fixed 
at the Galactic value ( 1 .06 x 10^" cm~^), estimated wit h the nH 
task of FTOOLS (based on lOickev & Lockman |[l990l) . 

We compute the effective area files (ARFs) and the re- 
sponse matrices (RMFs) for each region in the g rid. This pro- 




Fig. 11. Optical image with X-ray isocontours superimposed. 
The circles correspond to the three brightest galaxies in 
Abell 1758 North. The colors are the same as in Fig. \TQ\ the 
brightest galaxy is in black, the second brightest in blue and 
the third brightest in cyan. 



since we simultaneously fit all three instruments. The best fit 
value is then attributed to the central pixel. 

The resulting X-ray temperature and metallicity maps are 
displayed in Figs. [12] and [T3] Maps of the errors on these two 
maps are given in the Appendix (Fig. lA. 11 1. 

The temperature maps in Fig.[T2]show that both clusters ap- 
pear almost isothermal, the North one being hotter (~ 6-7 keV) 
and the South one being cooler (~ 4 - 5 keV). The presence of 
inhomogeneities is not obvious in the temperature maps, except 
for a hotter blob in the northwest region of Abell 1758 North. 
Although signs of recent interactions are not clearly present in 
the temperature map, the positions of the three brightest galax- 
ies argue in favour of a merging scenario where two clusters 
merged to form Abell 1758 North, as mentioned in Sect.|4] 

On the other hand, two striking aspects call the attention 
when considering the metallicity map shown in Fig. [13] First, 
we can note two elongated regions of high metallicity in the 
North cluster, suggesting that at least two s maller clusters have 



crossed the North cluster (as pointed out bv lDavid & Kempner 



cedure (already described in iDurret et al. I [201 01) allows us to 



perform a reliable spectral analysis in each spatial bin, in or- 
der to derive high precision temperature and metallicity maps. 



20041) . Based on the positions of the three brightest galaxies we 
can assume a scenario in which these elongated regions are due 
to metals were ton by ram-pressure stripping during the merg er 
Second, the metallicity map of the South cluster is even 
more unusual, since it shows a deficit of metals in the central 
region. This deficit is probably the signature of an interaction 
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Fig. 12. X-ray temperature map for Abell 1758 North (top) and 
South (bottom). The color-bar indicates the temperature in keV. 
The circles correspond to the three brightest galaxies in each 
cluster. The colors are the same as in Fig. [10] the brightest 
galaxy is in black, the second brightest in blue and the third 
brightest in cyan. The corresponding error map is shown in 
Fig.|Aj](top). 



33ml'^.4s 32m5^.f 



09.fis 13h31m48.0s 



ht! 

Fig. 13. X-ray metallicity map for Abell 1758 North (top) and 
South (bottom). The color-bar indicates the metallicity in solar 
units. The circles correspond to the three brightest galaxies in 
each cluster The colors are the same as in Fig. [10] the brightest 
galaxy is in black, the second brightest in blue and the third 
brightest in cyan. The corresponding error map is shown in 
Fig. lA.lK bottom). 



5. Numerical simulation results 



with the central object. Since this cluster is less massive, the ef- 
fects of galactic winds or supernova ex plosions will be stro nger 



lects or galactic wmds or supernova ex plosions will be stro 
to expell metals towards the outskirts (lEvrard et al. I2OO8I) 



However, none of the metal structures identified in the 2D 
maps correlates with any temperature structure for either clus- 
ter We do not see the high-metallicity regions of the North 
cluster in the temperature map. 

To have an overview of the merging scenario and better un- 
derstand the dynamical history of these clusters, we considered 
the results of numerical simulations to give support to our find- 
ings, and compared them with our observational results. This 
will be presented in the next section. 



David & Kempnerl ( 120041) have already analyzed this sys- 



tem using XMM-A^evvfo« and Chandra data. They conclude 
that these two clusters most likely form a gravitationally bound 
system, though their imaging and spectroscopic analyses do 
not reveal any sign of interaction between the North and South 
clusters, and show that A1758N and A1758S are both undergo- 
ing major mergers. It is clear from our analysis that Abell 1758 
North is still undergoing at least one merger, the bright and hot 
zone seen in the temperature map probably being the signature 
of one of the cores. For Abell 1758 South, the merger is not 
evident in X-rays, though some kind of interaction is suggested 
by the absence of a cool core, and confirmed by the optical data 
analysis. 



As des cri bed in detail in | Kapferer et al. I (l2006h iKapferer et al 



(l2007n or ISchindler et al. I (l2005b we combine several codes 
to model the metal distribution in a galaxy cluster To cal- 
culate the dark matter structure of the galaxy cluster we 



use the N-body code GADGET2 ( ISpringel 1 120051) with con 



strain ed random fields as initial conditions (IHoffman & Ribak 



1991]). A semi-analytical galaxy-formation model then assigns 



galaxy properties to the halos found in the dark matter 
structure. For this galaxy-formation mod el we use an im 



proved version of the code described by Ivan Kampen et al 



The chemical evolu ti on of the galaxies is modeled 



as in 



Matteucci & Fran9ois (Il989h . To study the evolution 



of the ICM we use a hydrodynamic grid code with comov- 
ing coo rdinates and a PPM-scheme for a better treatment of 
shocks ( IColella I I 1984. To optimize the computational time, 
we use 4 nested grids dRuffertll 19921) . each with 128x128x128 
cells. The total computational hydro-volume has a side length 
of 20 Mpc. We included special routines to transport en- 
riched material out from the modelled galaxies into the ICM: 
we model ram -pressure stripping as well as galactic winds 



as des c ribed in | Gunn & GottI (|1972|) : iDomainko & Kapferer 



(I2OO6I) : iKapferer et al. I (|2007|) . More specifically, the ram- 
pressure and galactic wind algorithms calculate a mass loss 
for each model galaxy (depending on its velocity relative to 
the ICM and its star formation rate). As the galaxy forma- 
tion model provides us with the metallicity within the ISM, 
we know the quantity of metals transported from the ISM into 
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the ICM. These metals are then advected with the ICM, and by 
knowing the ICM density we can then calculate the ICM metal- 
licity. Combining these codes, we can simulate the 3D evolu- 
tion of the ICM and its metallicity. From these 3D data, we 
then extract temperature and metal maps, by making emission 
weighted projections. These maps can be compared directly to 
observations. 

We performed 5 simulations with different initial condi- 
tions and selected the simulation which best matched our ob- 
servational results. It should be noted that we did not set up the 
simulations to specifically reproduce Abell 1758. Among all 
the tests perfomed here, we found one metal distribution that 
reproduces quite reasonably the elongated region of high abun- 
dance found observationally for the North cluster We see that 
the simulated 2D map in Fig. [14] can also reproduce the exact 
values for the metallicity of the North cluster For the bright 
regions (in yellow color) we have metallicity values of about 
0.5 - 0.6 solar units. The metallicity decreases towards the out- 
skirts where it reaches values around 0.2 solar units. 

It is of great importance to mention the spatial scales of the 
features presented in the 2D metal distribution maps predicted 
by numerical simulations and derived from observations. The 
physical sizes of the two top panels in Fig. [14] are exactly the 
same (a side-length of 2.5 Mpc), but the features, although sim- 
ilar in shape, do not have the same physical sizes. In the obser- 
vational metal distribution (top right panel in Fig. [T4l i we see 
that Abell 1758 North is completely comprised within 2.5 Mpc, 
while in the metal distribution computed from numerical sim- 
ulations (top left panel in Fig. [T4l) the cluster extends over a 
larger region. The point here is that we did not make the simu- 
lations to exactly match this cluster 

In order to analyze the importance of galactic winds and 
ram-pressure stripping in transporting metals, we display in 
Fig. [14] the results of numerical simulations, showing sepa- 
rately the metals ejected by galactic winds (bottom left), those 
transported out of the galaxies by ram-pressure stripping (bot- 
tom right), and the sum of the metals accounted for by these 
two processes (top left). These maps can be compared to the 
metallicity map derived from observations for the North cluster 
(top right). From these figures, we can say that both processes 
are important for metal enrichment, playing different roles in 
the cluster, with winds definitely playing a higher part in the 
south-west corner Galactic winds are more important in the 
outskirts of the cluster, while ram-pressure stripping seems to 
be dominant in the inner parts. 

It is worth noting the time scales of these different enrich- 
ment processes. By analyzing the results of our numerical sim- 
ulation, we can say that a large part of the metals for the high 
metallicity region in the south of Abell 1758 North and the thin 
elongated stripe from the center to the north-west have been 
transported into the ICM at redshifts z> 3. Then, between red- 
shifts 3 and 2 there was a substantial contribution to the metal- 
licity in the elongated structure but no contribution to the center 
of the cluster Between redshifts 2 and 1, there was a substantial 
contribution to the central region (around 1 /3 of the final metal- 
licity), mainly due to ram-pressure. Then, between redshifts 1 
and 0.8 there was hardly any contribution to the metallicity at 
all. From z=0.8 to z=0.5 there was a contribution (around 0.2 



solar metallicities) to the metals in the central region, coming 
only from ram-pressure stripping. And finally, from z=0.5 to 
z=0.279 (the cluster redshift) there was hardly any significant 
contribution to the metallicity. In summary, many of the metals 
which are responsible for the features in the outskirts have been 
transported into the ICM at very early times, even before z=3. 
If we look at the metals which have been transported into the 
ICM before z=2, we see that they explain the main morpho- 
logical features in the outskirts (like the elongated structure). 
On the other hand, the enrichment of the central region takes 
place mainly between z=2 and z=0.5, predominantly via ram- 
pressure stripping. 

The results presented above show that the comparison of 
observations with the results of numerical simulations is a pow- 
erful tool to understand better the physical processes involved 
in the transport of metals. However, there is still work to be 
done on this topic, based on a larger sample of clusters in dif- 
ferent dynamical states. 

6. Discussion and Conclusions 

It is becoming clear that clusters are in a continual state of 
evolution and with high precision telescopes, such as XMM- 
Newton and Chandra, we see that they are hardly ever relaxed 
and virialized as it was previously assumed. We have thus be- 
come interested in multiple mergers, where merging effects are 
expected to be even stronger 

As a second study of cluster pairs, we chose to analyse 
Abell 1758 (at a redshift of 0.279) by coupling archive optical 
CFHT Megacam and XMM-A^ewfo« data, and by comparing 
the temperature and metallicy maps obtained for the ICM with 
the results of numerical simulations. We made a step further by 
also considering the results of numerical simulations to try to 
understand better the dynamics and building-up history of this 
system. 

From our results, signs of merger(s) are detected in the op- 
tical as well as X-ray wavelength ranges, meaning that both 
galaxies and gas are still out of equilibrium. From optical re- 
sults we see that for Abell 1758 North, a Schechter function 
fits rather well most of the points of the GLFs in the g' band, 
but there is an excess of galaxies over a Schechter function 
in the brightest magnitude bins, specially in the r' band, as 
well as a possible and unexplained excess of galaxies around 
Mf ~ -17.5. On the other hand, for Abell 1758 South, which 
is poorer than the North cluster, the GLF is not as well defined 
and fit by a Schechter function. Note in particular that the GLFs 
derived by the two methods (statistical background subtraction 
for r' > 22 or colour-magnitude selection at all magnitudes)) 
look quite different, particularly in the g' band. If the dip de- 
tected in the r' band GLF for M,/ 18 is real, as already 

found in other clusters but at somewhat brighter absolute mag- 
nitudes, it could be interpreted as showing the tr ansition zone 



between (brighter) ellipticals and (fainter) dwarfs (iDurret et al 



An excess of galaxies over a Schechter function in the 
brightest magnitude bins is observed in both bands. We can also 
note that in the South cluster the GLFs in the g' and r' bands 
look quite different, with a faint end slope somewhat steeper in 
r' than in g'. The shallower faint end slope in g', if real, could 
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Fig. 14. Top left panel: metal distribution predicted by numerical simulations at the cluster redshift (i.e. sum of the two panels 
below). Top right panel: observed metal distribution map in Abell 1758 North. Both panels have a side-length of 2.5 Mpc. Bottom 
left panel: metals which were ejected by galactic winds. Bottom right panel: metals which were transported out of the galaxies 
by ram-pressure stripping. 



be explained either by quenching of star formation due to the 
merger (which strips galaxies from their gas and reduces star 
formation), or by the fact that star formati on has not ye t had 



time to be triggered by the merger (see e.g. iBekki 1119991 ). The 



somewhat perturbed shapes of the GLFs of both clusters there- 
fore agree with the fact that they are both undergoing merging 
processes. 

From the X-ray analysis we have noticed that the gas tem- 
perature maps do not present prominent inhomogeneities, ex- 
cept for a hotter blob in the northwest of Abell 1758 North. 
The North cluster is hotter (with temperatures in the range of 
~ 6 - 7 keV) and the South one cooler (with kT ~ 4 - 5 keV). 
The hotter blob in the northwest of Abell 1758 North could 
be explained by heating of the gas in that region by the move- 
ment of the northwest system towards the north proposed by 



David & Kempnerl (120041) . The most striking features are seen 
in the metallicity maps. The metallicity map of the South clus- 
ter is even more unusual, since it shows a deficit of metals in the 
central region. This deficit is probably the signature of an inter- 
action with the central o bject that could hav e expelled metals 
towards the outskirts (see lEvrard et al. 1l2008[ and top left panel 
of Fig.[T4li. We also detect two elongated regions of high metal- 
licity in the North cluster, suggesting that at least two smaller 
clusters have crossed the North cluster Note that our tempera- 
ture and met alhcity maps are in agreement with the scenarios 
proposed by lOavid & Kempner I (|2004), who derived that the 
North cluster is in the later stages of a large impact parameter 
merger between two 7 keV clusters, while the south cluster is 
in the earlier stages of a nearly head-on merger between two 
5 keV clusters. 
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In order to understand better the nature of the most promi- 
nent features exhibited in the metallicity map of the North clus- 
ter, we performed 5 simulations with different initial condi- 
tions. It should be stressed that we did not set up the simu- 
lations to specifically reproduce Abell 1758. Among these 5 
simulations, we found one metal distribution that reproduces 
quite reasonably the elongated region of high abundance found 
observationally for the North cluster, although without a per- 
fect spatial correlation. The results of our numerical simula- 
tions allowed us to distinguish the role of metal transportation 
processes such as galactic winds and ram-pressure stripping. 
We have shown that these phenomena act in different regions 
of the cluster, and that in the metal-rich elongated regions of 
the North cluster, winds were more efficient in transporting en- 
hanced gas to the outskirts than ram pressure stripping. These 
simulations also allowed us to compute the GLF, and the result 
is consistent with the observed GLFs. 

In the hope of finding some kind of large scale structure 
and/or filaments linking Abell 1758 with its surroundings, we 
searched the NED database for galaxies with redshifts available 
in a region of 36 arcmin around Abell 1758 (corresponding to 
9.1 Mpc at the cluster distance). We found 137 galaxies, out of 
which only 15 have redshifts in the approximate cluster redshift 
range: 0.264, 0.294 . No structure is found in the blueshifted 
and redshifted galaxies, except for what appears to be a cluster 
west of Abell 1758 (see Section 3.2). Therefore we cannot de- 
rive any argument from the large scale distribution of galaxies 
around Abell 1758. 

To conclude, the Abell 1758 North and South clusters most 
likely form a gravitationally bound system, but our imaging 
analyses of the X-ray and optical data do not reveal any sign 
of interaction between the two clusters. All signs of dynamical 
disturbance are associated with recent merger(s) which each of 
the two clusters is still undergoing. 
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Appendix A: Error maps 

In Eig. lA.ll we present the errors associated to each bin in the 
temperature and metallicity maps displayed in Eigs.fTSjandfTSl 
The errors on these parameters were directly obtained from the 
spectral fits. 

Looking at Ei g s . [T2l 1 1 3 1 and IATTI (error maps) we see that, 
for the temperature estimates, we have errors around 5% for 
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Metallicity error map 
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Fig.A.l. Upper panel: error map on the X-ray temperature. 
Lower panel: error map on the X-ray metallicity. 



A 175 8 South, while for A 175 8 North they vary from 5% (in 
the inner parts) up to 12% (at the outskirts). The errors on the 
temperature maps are therefore very reasonable. For the metal- 
licity maps we do not have the same accuracy. Looking at the 
same figures, we see that for A1758 South, the metallicity er- 
rors can reach 25% and for A1758 North they vary from 30% 
(in the inner region) up to 46% (at the outskirts). It is impor- 
tant to notice that although the XMM-Newton exposure time 
on Abell 1758 was about 57 ks, due to flare filtering only a 
17 ks of good data were useable to construct the 2D spectral 
maps. 



